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Abstract —In this paper, we introduce the novel use of linear 
spatial precoding based on fixed and known parameters of 
multiple-input multiple-output (MIMO) channels to improve the 
performance of space-time coded MIMO systems. We derive 
linear spatial precoding schemes for both coherent (channel 
is known at the receiver) and non-coherent (channel is un¬ 
known at the receiver) space-time coded MIMO systems. Antenna 
spacing and antenna placement (geometry) are considered as 
fixed parameters of MIMO channels, which are readily known 
at the transmitter. These precoding schemes exploit the antenna 
placement information at both ends of the MIMO channel to 
ameliorate the effect of non-ideal antenna placement on the 
performance of space-time coded systems. In these schemes, 
the precoder is fixed for given transmit and receive antenna 
configurations and transmitter does not require any feedback 
of channel state information (partial or full) from the receiver. 
Closed form solutions for both precoding schemes are presented 
for systems with up to three receiver antennas. A generalized 
method is proposed for more than three receiver antennas. We 
use the coherent space-time block codes (STBC) and differential 
space-time block codes to analyze the performance of proposed 
precoding schemes. Simulation results show that at low SNRs, 
both precoders give significant performance improvement over a 
non-precoded system for small antenna aperture sizes. 

Index Terms —Space-time coding, channel modelling, linear 
precoder design, MIMO systems, non-isotropic scattering, spatial 
correlation. 


I. Introduction 

M IMO communication systems that use multi-antenna ar¬ 
rays simultaneously during transmission and reception 
have generated significant interest in recent years. Under the 
assumption of fading channel coefficients between different 
antenna elements are statistically independent and fully known 
at the receiver (coherent detection), theoretical work of [1] 
and [2] revealed that the channel capacity of multiple-antenna 
array communication systems scales linearly with the smaller 
of the number of transmit and receive antennas. Motivated 
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by these works, [3-5] have proposed several modulation and 
coding schemes, namely space-time trellis codes and space- 
time block codes, to exploit the potential increase in capacity 
and diversity gains using multi antenna arrays with coherent 
detection. 

The effectiveness of these coherent space-time coding 
schemes mainly relies on the accuracy of the channel es¬ 
timation at the receiver. Therefore, differential space-time 
coding (DSTC) schemes proposed in [6-8] make an attractive 
alternative to combat inaccuracy of channel estimation in 
coherent space-time coding schemes. With DSTC schemes, 
channel state information is not required at either end of 
the channel. However, it is well known that DSTC schemes 
suffer a 3dB performance loss compared to space-time coding 
schemes with coherent detection at the receiver. 

For both schemes, code structures are designed assuming 
that the channel gains between the transmitter and receiver 
antennas undergo uncorrelated independent flat fading. Such 
an assumption is valid only if the scattering environment 
is isotropic, i.e., scattering is uniformly distributed over the 
receiver and transmitter antenna arrays, and also only if the 
antennas in an array are well separated. Recent studies have 
shown that insufficient antenna spacing and non-isotropic scat¬ 
tering reduce the performance of space-time coded communi¬ 
cation systems [9-11]. This has motivated the design of linear 
precoders for space-time coded multiple antenna systems by 
exploiting the statistical information of the MIMO channels 
[12-17]. In these schemes, the receiver either feeds back the 
full channel state information or the correlation coefficients 
of the channel (covariance feedback) to the transmitter via 
a low rate feedback channel. In order to be cost effective 
and optimal, these designs assumed that the channel remains 
stationary (channel statistics are invariant) for a large number 
of symbol periods and the transmitter is capable of acquiring 
robust channel state information. However, when the channel 
is non-stationary or it is stationary for a small number of 
symbol periods, the receiver will have to feedback the channel 
information to the transmitter frequently. As a result, the 
system becomes costly and the optimum precoder design, 
based on the previously possessed information, becomes out¬ 
dated quickly. In some circumstances feeding back channel 
information is not possible. These facts have motivated us 
to design a precoding scheme based on fixed and known 
parameters of the underlying MIMO channel. Following list 
summarizes the original contributions of this paper. 

• By exploiting the spatial dimension of a MIMO channel. 


we design linear spatial precoding schemes to improve 
the performance of coherent and differential space-time 
block coded systems. These linear spatial precoders are 
designed based on previously unutilized fixed and known 
parameters of MIMO channels, the antenna spacing and 
antenna placement details. We use the spatial channel 
decomposition given in [18] to incorporate the antenna 
spacing and antenna placement details into the precoder 
design. 

• Both precoders are fixed for fixed 1 antenna placement and 
the transmitter does not require any form of feedback 
of channel state information (partial or full) from the 
receiver. 

• Since the designs are based on fixed parameters, these 
spatial precoders can be used in non-stationary channels 
as well as stationary channels. 

• Upper bounds for pairwise error probability (PEP) of 
coherent space-time codes and differential space-time 
codes are derived for spatially correlated MIMO chan¬ 
nels. To the authors knowledge, the PEP upper bound 
of differential space-time codes is a new bound. Uti¬ 
lizing the MIMO channel decomposition given in [18], 
antenna configuration details and scattering environment 
parameters (angular spreads and mean angle of arrival 
and departure) are incorporated in to these PEP upper 
bounds. Assuming an isotropic scattering environment 
surrounding the transmitter and receiver antenna arrays, 
we minimize the two PEP upper bounds to obtain the 
optimum precoders. 

• We show that our spatial precoding schemes reduce the 
effect of non-ideal antenna placement, which is a major 
contributor to the spatial correlation, on the MIMO sys¬ 
tem performance. In these schemes, the precoder virtually 
arranges the antennas into an optimal configuration as 
such the spatial correlation between all antenna elements 
is reduced. 

• The precoder design is based on the spatial channel model 
proposed in [18], but we show that the performance of 
both precoding schemes does not depend on the channel 
model that used to model the underlying MIMO channel. 
Therefore, our design and simulation results provide an 
independent confirmation of the validity and usefulness 
of the channel model proposed in [18]. 

An outline of the paper is as follows. Section previews the 
spatial channel model used in our design. In Section |TO] the 
precoded coherent STBC and differential STBC systems are 
described along with detection rules at the receiver. Sections 
EYI and El present the optimization problem and the optimal 
precoder solution for coherent STBC and differential STBC, 
respectively. For both precoding schemes, we show that the 
optimum linear precoder for a multiple-input single-output 
(MISO) fading channel is essentially given by the classical 
“water-filling” strategy found in information theory [1]. For a 
MIMO channel, the linear precoder is determined by a novel 
generalized water-filling scheme. Closed form solutions for 
both precoding schemes are presented for systems with up 

'antennas are fixed relative to each other 


to three receiver antennas. A generalized method is proposed 
for more than three receiver antennas. Sections eh and Em 
present results obtained with proposed precoding schemes 
for various spatial scenarios using the spatial channel model 
in [18] as the underlying MIMO channel. Section EIH also 
presents results obtained with proposed precoding scheme 
for non-isotropic scattering environments (i.e., limited angular 
spreads at the transmitter and receiver antenna arrays). Section 
IvTm gives the simulation results of our proposed precoding 
scheme applied on other statistical channel models found in 
the literature. Section llxl present some concluding remarks and 
five appendices contain various proofs. 

Notations: Throughout the paper, the following notations 
will be used: Bold lower (upper) letters denote vectors 
(matrices). [-] T , [•]* and [■]' denote the transpose, complex 
conjugate and conjugate transpose operations, respectively. 
The symbols S(-) and <g> denote the Dirac delta function and 
Matrix Rronecker product, respectively. The notation E {■} 
denotes the mathematical expectation, vec(A) denotes the 
vectorization operator which stacks the columns of A, tr{ -} 
denotes the matrix trace, [.] denotes the ceiling operator and 
S 1 denotes the unit circle. The matrix I n is the nxn identity 
matrix. 


II. Spatial Channel Model 
First we review the spatial channel model proposed in 
[18]. Consider a MIMO system consisting of nr transmit 
antennas located at positions u t , t = 1, 2, • • • ,ut relative to 
the transmitter array origin, and hr receive antennas located 
at positions v r , r = 1,2, ••• ,ur relative to the receiver 
array origin, tt > max || u t || and tr > max || v r || 
denote the radius of spheres that contain all the transmitter 
and receiver antennas, respectively. We assume that scatterers 
are distributed in the far field from the transmitter and receiver 
antennas and regions containing the transmit and receive 
antennas are distinct. 

By taking into account physical aspects of scattering, the 
MIMO channel matrix H can be decomposed into determin¬ 
istic and random parts as [18] 

H = JrHsJU (1) 


where Jr is the deterministic receiver configuration matrix, 


J—Nr (^ i ) 
J-Nr(v 2 ) 


R (til ) 

Jn r (v 2 ) 


3— Nr { v ur ) JNr {^riR ) 


and Jr is the deterministic transmitter configuration matrix. 


J-n t ( u 1 ) 
J-n t ( u 2 ) 


Jn t ( u 1 ) 
Jn t { u 2) 


J—Nt i^ , n r p ) 


Jnt irp ) 


J ri (w ) is the spatial-to-mode function (SMF) which maps the 
antenna location w to the n-th mode of the region. The form 






which the SMF takes is related to the shape of the scatterer- 
free antenna region. For a circular region in 2-dimensional 
space, the SMF is given by a Bessel function of the first kind 
[18] and for a spherical region in 3-dimensional space, the 
SMF is given by a spherical Bessel function [19]. For a prism¬ 
shaped region in 3-dimensional space, the SMF is given by a 
prolate spheroidal function [20]. 

Here we consider the situation where the multipath is 
restricted to the azimuth plane only (2-D scattering envi¬ 
ronment), having no field components arriving at significant 
elevations. In this case, the SMF is given by 

where J n (■) is the Bessel function of integer order n, vector 
w = (|| w ||, (j> w ) in polar coordinates is the antenna location 
relative to the origin of the aperture, k = 2tt/X is the wave 
number with A being the wave length and i = \/ - 1. Jr is 
riT x (2Nr -f 1) and Jr is urx(2Nr + 1), where 2Nt + 1 and 
2 Nr + 1 are the number of effective 2 communication modes 
at the transmit and receive regions, respectively. Note, Nr and 
Nr are defined by the size of the regions containing all the 
transmit and receive antennas, respectively [21]. In our case, 

and 


where e « 2.7183. 

Finally, H$ is the (2Nr + 1) x (2Nr + 1) random complex 
scattering channel matrix with th element given by 

{H s } e m = ff g(4>, 

J J^x S 1 

( 2 ) 

representing the complex scattering gain between the (m — 
Nt — l)-th mode of the scatter-free transmit region and 
(£ — Nr — l)-th mode of the scatter-free receiver region, 
where g{(f >, ip) is the effective random complex scattering gain 
function for signals with angle-of-departure (f> from the scatter- 
free transmitter region and angle-of-arrival <p at the scatter-free 
receiver region. 

The channel matrix decomposition 0 separates the chan¬ 
nel into three distinct regions of interest: the scatter-free 
region around the transmitter antenna array, the scatter-free 
region around the receiver antenna array and the complex 
random scattering environment which is the complement of 
the union of two antenna array regions. Consequently, the 
MIMO channel is decomposed into deterministic and random 
matrices, where deterministic portions Jr and Jr represent 
the physical configuration of the transmitter and the receiver 
antenna arrays, respectively, and the random portion represents 
the complex scattering environment between the transmitter 
and the receiver antenna regions. The reader is referred to 
[18] for more information regarding this spatial channel model. 


Nr — 
Nr = 


kerr 

A 

ker R 
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Note that the precoder design is based on this channel model, 
but the performance does not depend on this model (see 
Section ivmi . That is, our design and simulations provide an 
independent confirmation of the validity and usefulness of this 
channel model. 

A. Spatial Correlation 

Suppose transmitter configuration matrix Jr has the singu¬ 
lar value decomposition (svd) Jr = UrArV' t and receiver 
configuration matrix Jr has the svd Jr = U rArV r . 
Substituting svds of Jr and J r in 0 and using the Kro- 
necker product identity [22, page 180] vec {AXB) = ( B T ® 
A) vec(X), we obtain 

h = hjs(U T R ®U^ T ), (3) 

where hjs = (vec(iTj S )) T with Hjs = 
ArV^HsVrA^. Applying the same Kronecker product 
identity to vec (H T JS ) yields hj$ = ^s , [(V'^AU)®>(VtA^)], 
where hs = (vec(iTj)) 7 . Then the covariance matrix Rjj 
of the MIMO channel H is given by 

Rh = £{h + h}, 

= {U*r®U t )Rjs{U T r®U^ t ), (4) 

where R JS = [(A^Vj) ® (A t V^ t )]R s [{V* r aI) <g> 

(VrA^)} with R s = E {h^h 5 }. 

In this work, our main objective is to design a linear 
precoder which compensates for any detrimental effects of 
non-ideal antenna placement/configuration on the performance 
of space-time block codes. Here we assume that the scatter¬ 
ing environment surrounding the transmitter and the receiver 
regions is “rich 3 ”, i.e., R§ = I. This assumption yields the 
simplification 

Rjs = [(A * r V T r) ® (A T V^)][(n.A£) ® (VrAt,)] (5a) 
= (A* R Al) ® (A t a£), (5b) 

where m follows from (l5al by matrix identity [22, page 180] 
(. A ® C)(B ® D) = AB ® CD, provided that the matrix 
products AB and CD exist, and unitary matrix properties 
V r Vr = I and V r Vr = I. Substituting J5bl into 0 gives 

Rh = {U* r ®Ut){Rr®Rt)(U t r ®U^ t ), (6) 

where 

Rr = ArA jp (7) 

and 

Rr = (ArA^ r ) t . (8) 

Note that both Rr and Rt are diagonal matrices, where the 
diagonal of Rr consists of squared singular values of Jr 
(or eigen-values of JrJ^ r ) and diagonal of Rt consists of 
squared singular values of Jr (or eigen-values of JrJr)- 


2 Although there are infinite number of modes excited by an antenna array, 
there are only finite number of modes (27V + 1) which have sufficient power 
to carry information. 


3 Even though precoders are derived for rich scattering channels, these 
precoders provide significant perfo rmance improvements in non rich scattering 
channel environments, see Section Ivii-cl 






III. System Model 

At time instance k, the space time encoder at the 
transmitter takes a set of modulated symbols C(k) = 
{ci(fc), C 2 (fc), • • • ,Cft-(/c)} and maps them onto an nrxT 
code word matrix Sgr k ) 6 V of space-time modulated con¬ 
stellation matrices set V = {Si, S 2 , • • • , Sl}, where T is the 
code length, L = q K and q is the size of the constellation 
from which c n (k), n = 1, • • • ,K are drawn. By setting 
\c n (k)\ = 1/y/K, each code word matrix S^ k j in V will 
satisfy the property S^ k )S' e , k , = I nT for £(k) = 1,2, ■ • • , L. 

In this paper, we mainly focus on the space-time modulated 
constellations with the property 

(Si - Sj)(Si - srf = 0iji nT yi ? j, (9) 

where is a scalar and Si, Sj £ V. Space-time orthogonal 
designs in [5] and some cyclic and dicyclic space-time mod¬ 
ulated constellations in [7] are some examples which satisfy 
property © above. 

A. Coherent Space-time Block Codes 

Let s n be the n-th column of Si = [si, S 2 , • • • , St] £ V- At 
the transmitter, each code vector s n is multiplied by a njx 
nr fixed linear precoder matrix F c before transmitting out 
from nr antennas. Assuming quasi-static fading, the signals 
received at tir receiver antennas during T symbol periods can 
be expressed in matrix form as 

Y(k) = s/W s HF c S l[k) + N(k ), 


with X(0) = I nT . Then, each encoded X(k) is multiplied by 
a riT^nr fixed linear precoder matrix F d before transmitting 
out from nr transmit antennas. Assuming quasi-static fading, 
the signals received at hr receiver antennas during nr symbol 
periods can be expressed in matrix form as 

Y(k) = s/E~ s HF d X(k) + N(k), 

where N(k) is the riRXriT white Gaussian noise matrix 
in which elements are zero-mean independent Gaussian dis¬ 
tributed random variables with variance a 2 /2 per complex 
dimension and H is the hr x tit channel matrix, which is 
modeled using 0 . 

Assume that the scattering channel matrix Hr remains 
constant during the reception of two consecutive received 
signal blocks Y(k — 1) and Y(k), which can be expressed 
in vector (row) form as 

y(k — 1) = \/~E(hX(k — 1) + n(k — 1), 
y(k) = s/WshX (k) +n(k), 

= y(k~l)S m +w(k), (11) 

where y(k ) = (vec (Y(k) T )) , X(k) = (g> (F d X(k)), 

h = (vec (H T )) T , n(k) = (vec (N(k) T )) , S e ( k ) = In R ® 
S e (k) and w(k ) = n(k ) - n(k - 1 )S m . 

For differential STBC, we assume that receiver has no CSI 
whilst transmitter has partial CSI (antenna placement details). 
From the transmitted code word matrix is detected 

differentially using the minimum Euclidian distance detection 
rule: 


where E s is the average transmitted signal energy per symbol 
period, N(k) is the urxT white Gaussian noise matrix 
in which elements are zero-mean independent Gaussian dis¬ 
tributed random variables with variance cr ^/2 per dimension 
and H is the n r x ut channel matrix. In this work, we use the 
channel decomposition 0 to represent the underlying MIMO 
channel and the elements of scattering channel matrix Hr 
are modeled as zero-mean complex Gaussian random variables 
(Rayleigh fading). 

For coherent STBC, we assume that the receiver has perfect 
channel state information (CSI) and transmitter has partial CSI 
(antenna placement details). At the receiver, the transmitted 
codeword is detected by applying the minimum Euclidian 
distance detection rule: 

St(k) = arg min || y(k) - \/E~ s hS^ k ) || 2 
Se(k)€V 

= arg max Rc{hS^ k ) y f (fc)}, (10) 

S , £(i s )SV 

where y(k) = (vec (Y T (k))) T , S^ k ) = I nR <g) S t ( k) and 
h = (vec (H T )) with H = HF C . 

B. Differential Space-time Block Codes 

In this scheme, codeword matrix Sg( k ) is differentially 
encoded according to the rule 

X(k) = X(k-l)S m , for fc = 1, 2, • • • 


Se(k) = arg 


min 

^(fc)GV 


y(k) - y(k - 1 )S e{k) || 2 


= arg max R e{y(k - l)S m y(k ) 1 }. 
■S j(t)6V 


IV. Problem Setup: Coherent STBC 


Assume that perfect CSI is available at the receiver and 
also maximum likelihood (ML) detection is employed at the 
receiver. Suppose codeword Si £ V is transmitted, but the 
ML-decoder m chooses codeword Sj £ V, then as shown 
in the AppcndixQ] the average pairwise error probability (PEP) 
is upper bounded by 


P(Si - Sj) < 


IRnlln 


s a] 


( 12 ) 


where S A = F c (Si - Sj)(S t - Sj^Fl, R h = E |h f h} 

with row vector h = (vec (H t )) t and 7 = E s /cr R is the 
average symbol energy-to-noise ratio (SNR) at each receiver 
antenna. Substituting 0 in (1 1 21 and applying the property 0 
associated with orthogonal space-time block codes we obtain 


P (Si - Sj) < 


1 

In T n R + R RT [In R 0 U\,F c F\U T } ’ 

(13) 


where we have used the matrix determinant identity 
\I + AB\ = 1 1 + BA\ and let Rrt = Rr 0 Rt 









Optimization Problem 1: Find the optimum spatial precoder 
F c that minimizes the average PEP upper bound (1131 for 
coherent STBC, subject to the transmit power constraint 
trl-Fc-F^} = tit , for given transmitter and receiver antenna 
configurations in a rich scattering environment. 


A. Optimum Spatial Precoder: Coherent STBC 

The linear precoder F c is designed by minimizing the 
maximum of all PEP upper bounds subject to the power 
constraint tr{.F c .F c '} = tit - Alternatively, let 


Q c= lJh±U^ r F c F\U T , 

then the average PEP bound O becomes 

1 


P (Si - Sj) < — 

J- n 


+ [R R ®R T ][In R ®Q c }\ ' 


(14) 


and Q c must satisfy the power constraint tr{ Q r } = 

^T7/3fc,^/4. Since log(-) is a monotonically increasing func¬ 
tion, the logarithm of the average PEP upper bound O is 
used as the objective function to minimize. Note that Q c in 
ill41 is always positive semi-definite as Q c = BB\ with 
B = V(7/3m)/4 U\F C . 

Now the optimum Q c is obtained by solving the optimiza¬ 
tion problem: 


min - log \In T n R + (Rr ® Rr)(In R ® Qc)\ 

subject to Q c R 0, tr{Q c } = (15) 

where j3 = min k^l{Pk,t} over a ll possible codewords 4 . By ap¬ 
plying Hadamard’s inequality on \I + (Rr ® Rt)(I ® Q c )\ 
gives that this determinant is maximized when (Rr®Rt)(I® 
Q c ) is diagonal [1]. Therefore Q c must be diagonal as Rr and 
Rt are both diagonal. Since (Rr® Rt)(I ®Q c ) I s a positive 
semi-definite diagonal matrix with non-negative entries on its 
diagonal, I + (Rr ® Rt)(I ® Q c ) forms a positive definite 
matrix. As a result, the objective function of our optimization 
problem is convex [23, page 73]. Therefore the optimization 
problem 03 above is a convex minimization problem because 
the objective function and inequality constraints are convex 
and equality constraint is affine. 

Let qi = [Q c ]i,i , U = [R T \i,i and rj = [RrIjj- Optimiza¬ 
tion problem 03 then reduces to finding qi> 0 such that 


min — EE log(l + tiqiVj) 

j=i i =i 

subject to q R 0, 

1 T g = n TlP 


(16) 


where q = [< 71 , < 72 , • • • j qn T ] T and 1 denotes the vector of all 
ones. 

Introducing Lagrange multipliers A c S R'" 7 for the inequal¬ 
ity constraints — q R 0 and v c 7 ]R. for the equality constraint 


l 7 q = tit"/3 /4, we obtain the Karush-Kuhn-Tucker (K.K.T) 
conditions 


q h 0, A c R 0, 1 7 q = 

Kqt — 0; l — I? 2, • 


n.T'yfi 

4 

, n T 


-St 

i =1 


r-iti 


+ rjUq 


Aj T v c — 0, i — 1,2, 


,n T . (17) 


A i in d 1 7> can be eliminated since it acts as a slack variable 5 , 
giving new K.K.T conditions 


qh 0 , i J <7 = 


riT7/3 


Qi Vc 


-E 

3 =1 
riR 

^E 

3 =1 


r iU 


1 + rjUqi 

r jtj 

1 + rjUqi 


= 0, * = n T , 


, i = !,■•• , n T - 


(18a) 

(18b) 


For hr = 1, the optimal solution to 03 is given by the 
classical “water-filling” solution found in information theory 
[1]. The optimal q, for this case is given in Section II V- 131 
For hr > 1, the main problem in finding the optimal qi for 
given ti and rq,j = 1,2,- •• ,ur is the case that, there are 
multiple terms that involve ^ on fTO . Therefore we can view 
our optimization problem 03 as a generalized water-filling 
problem. In fact the optimum q.j for this optimization problem 
is given by the solution to a polynomial obtained from ill Sat . In 
Sections EH] and IIV-DI we provide closed form expressions 
for optimum qi for hr = 2 and 3 receiver antennas and a 
generalized method which gives optimum qi for hr > 3 is 
discussed in Section HV-EI 

As shown above, the optimal Q c is diagonal with 


Q c = diag{<7i,<72, ■ • • ,q nT }, 


and optimal spatial precoder F c is obtained by forming 



where U n is any unitary matrix. In this work, we set U n = 

ITIT ’ 


B. MISO Channel 


Consider a MISO channel where we have riT transmit an¬ 
tennas and a single receive antenna. The optimization problem 
involved in this case is similar to the water-filling problem in 
information theory, which has the optimal solution 

qi = { Vc<ti ’ 

l 0, otherwise, 

where the water-level l/v c is chosen to satisfy 



4 


4 Setting (3 = minwill minimize the error probability of the 
dominant error event(s). 


5 If g(x) < v is a constraint inequality, then a variable A with the property 
that g(x) + A = v is called a slack variable [23]. 




















V. Problem Setup: Differential STBC 


C. riT X 2 MIMO Channel 


We now consider the case of ut transmit antennas and 
na = 2 receive antennas. As shown in the Appendix [H] the 
optimum q, for this case is 


Qi = 


A + VK, v c < U(ri + r 2 ); 

0 , otherwise. 


(19) 


where v c is chosen to satisfy 


max (o, A + Vl?) = 


with 


A = 


K = 


2r 1 r 2 t 2 i - VcUiri +r 2 ) 
2v c r 3 r 2 t 2 

\t 2 (r i - r 2 ) 2 + 4 r\r 2 tj 
2v c n r 2 t 2 


and 


( 20 ) 


D. riT x 3 MIMO Channel 

For the case of tit transmit antennas and n r = 3 receive 
antennas, the optimum qi is given by 


/ ~S + ‘ S ' + T, v c < ti(n + r2 + r 3 ); 

\ 0 , otherwise. 


( 21 ) 


where v c is chosen to satisfy 



0, + S + T 

Sa 3 


riT'yP 

4 


with 

S + T = 


R + VQ 3 + R 2 


^ _ 3 aia 3 — a 2 
9a| ’ 


R = 


R - VQ 3 + R 2 \ , 
9010203 — 27aga 2 — 2 ai\ 


54 a\ 


03 = v c rir 2 r 3 ti, a 2 = v c t 2 (rir 2 + 77 r 3 + r 2 r 3 ) - 37Tr 2 r 3 <f, 
ai = v c ti(r i + r 2 + r 3 ) - 2 t 2 (nr 2 + r ± r 3 + r 2 r 3 ) and 
ao = v c — ti(ri + r 2 + r 3 ). A sketch of the proof of (1211 is 
given in the Anpendix-Hnl 


E. A Generalized Method 

We now discuss a method which allows to find optimum 
solution to for a system with nr transmit and hr receive 
antennas. The complementary slackness condition X,q, = 0 
for i = 1,2, ••• ,tit states that A * is zero unless the i-th 
inequality constraint is active at the optimum. Thus, from 
(II 8 al we have two cases: (i) = 0 for v c > t,; i Tj> 

(ii) v c = r j i </(l + r .-jUqi) for qi > 0 [23, page 243]. 
For the later case, the optimum qi is found by evaluating the 
roots of hr -th order polynomial in q ,, where the polynomial 
is obtained from v c = Y^ij=i Tjt*/(1 + r jtiQi). Since the 
objective function of the optimization problem is convex 
for q > 0, there exist at least one positive root to the riR-th 
order polynomial for v c < U i r r l fl the case °f multiple 
positive roots, the optimum qi is the one which gives the 
minimum to the objective function of (1161 . In both cases, v c 
is chosen to satisfy the power constraint 1 1 q = n-r 7 / 3 / 4 . 


For the Differential STBC, we again use the average PEP 
upper bound to derive the optimum spatial precoder that 
reduces the effects of non-ideal antenna placement on the 
performance of differential STBC. Below shows the derivation 
of the average PEP upper bound. 

Based on (DU, the receiver will erroneously select Sj when 
Si was actually sent as the fc-th information matrix if 

II y(k) - y{k - 1 )Sj || 2 < || y(k) - y{k - 1)5, || 2 , 

y{k- l)Dijy\k- 1) < 2Re{w(k)A\^ (k - 1 )}, 

( 22 ) 


where A,;.,- = Sj — S, = I nR <g> (S :j - S t ) and D it j = 
Aij At . = I^dSi-Sj^Si - SjY). For given y(fc-l), 
the term on the left hand side of <E3 is a constant and the 
term on the right hand side is a Gaussian random variable. Let 
u = 2Re{ut(fc)AT ^{k — 1)}, then in the Appendix IIVI we 
have shown that u has the conditional mean 


m u \y(k-i) = E{u\y{k- 1)} , 

= 2Re |m T i(fc-i)| J/ (fc-i)(T - S l S ] j )y\k - 1)|, 

where = a 2 n y{k - 1 ){X\k - 1 )R H X(k - 

1 ) + o' 2 J„ Tnjl ) _1 , and the conditional variance 

a t\ y (k- 1 ) = E {\\u- rh u \ y (k-i) || 2 | y{k - 1)} , 

= 2 y(fc - 1 )A itj 

X nI T 1)|l/(fc— 1)^*) ~ 1)) 

where E„( fe _i)| J/(fe _i) = tj 2 n (I - cj 2 n (E s X\k - l)R H X(k - 
1) + ct 2 J) _1 ). Recall that R H in rh n (k-i)\ y (k-i) and 
’En(k-i)\y(k-i) the channel correlation matrix, defined by 
0 and X{k) = I nR ® (F d X(k)). 


Let d 2 j = y(k — 1 )Dijy^(k — 1). Based on (1221 . the PEP 
condition on received signal y(k — 1 ) is given by 


P{Si -> Sj | y{k - 1 )) 


Pr(t/ > dlj), 


d? 


\f2tta 


■ exp 



(u — to ) 2 

2 <t 2 


du, 


(23) 


In order to obtain unconditional PEP, we need to average o 
with respect to the distribution of y(k — 1). Unlike in the 
coherent STBC case, finding unconditional PEP from o 
poses a much harder problem due to the non-zero fh u \ y (k-i) 
and complicated cr 2 | i/ ^ fc _ 1 ^. However, at asymptotically high 
SNRs (i.e., keeping E s constant and cr 2 ^ 0) the conditional 
mean and the conditional variance of u reduce to TO u | 1 / (fc_i )=0 
and = 4 a R d 2 j, respectively. As shown in the 

Appendix [V] the average PEP can be upper bounded by 


P(S( - Sj) < 

1 


/ + i (jX(k - 1)^R h X{ k - 1) + I 


nrriR I D ij 


( 24 ) 




















where 7 = E s /a 2 is the average SNR at each receiver 
antenna. As for the coherent STBC case, we mainly focus 
on the space-time modulated constellations with the property 
(|9j. Furthermore, similar to [7, 8 ] we assume that code length 
T = nr . Under this assumption, each code word matrix £7 in 
V will satisfy the unitary property S,S] = I and S ' t S', = I 
for i = 1 , 2, • • • , L. As a result, X(k) will also satisfy the 
unitary property X(k)X^ {k) = I and x\k)X(k) = I for 
k = 0 , 1,2, •••. Applying 0 on and then using the 
unitary property of X(k — 1 ) and the determinant identity 
1 1 + AB | = \I + BA |, after straight forward manipulations, 
we can simplify the PEP upper bound (1241) to 


P(S. - S s ) < - 


(s+|ii) 


( 8 +0i,j) 


Rh^ur <S> F d Fj) 


(25) 


As before, we assume that the scattering environment 
surrounding the transmitter and receiver antenna arrays is 
isotropic. Then, substitution of © in gives 


P(Si - Sj) < 


1 

2 


8 + gi,.A-" T "* 

8 ) _ 

I + tM 7 )( Rr ® ® U' T F d F\UT) 


(26) 


where Rt and Rr are defined by 0 and 0 , respectively. 
The optimization problem for differential STBC case can 
now be stated as follows: 


Optimization Problem 2: Find the optimum spatial precoder 
F d that minimizes the average PEP upper bound d26l for 
differential STBC, subject to the transmit power constraint 
ti{F d F^ d } = riT, for given transmitter and receiver antenna 
configurations in a rich scattering environment. 


A. Optimum Spatial Precoder: Differential STBC 

Similar to the coherent STBC case, the optimum spatial 
precoder F d for differential STBC is obtained by minimizing 
the maximum of all PEP upper bounds subject to the power 
constraint trjF’dF 1 ^} = tit- Let 


Pd = 


Pi,o 7 

(8 + fiij) 


U ] T F d F\u T . 


The optimum P,i (hence the optimum Ff) is then obtained 
by solving the optimization problem 


min -log 11 + (Rr ® R T ){I nR 0 Pd)\ 
subject to P d A 0, tr{P d } = ■ 

(8 + p it j) 

The above optimization problem is identical to the optimiza¬ 
tion problem derived for coherent STBC, except a different 
scalar for the equality constraint. Therefore, following Section 
EMI here we present the final optimization problem and 
solutions to it without detail derivations. 


Following Section llV-AI we can show that the optimum P,i 
is diagonal and diagonal entries of P,i are found by solving 
the optimization problem 


riR n T 

min — EE log(l -) -tiPiCj) 

j=l i =1 


subject to p A 0 , 


1 T P 


fPyriT 

(8 + /3) 


(27) 


where /3 = min,;^ 7 {/3,; >; } over all possible codewords. 
Pi = \Pd]i,i, U = [Rr]i,i rj = [RrIjj and p = 
[pi,P 2 , ■ ■ ■ ,Pn T ] T - The linear spatial precoder F d is obtained 
by forming 


F d = 



where P d = diag{pi,p 2 , ■ ■ • ,Pn T } an d Ff n is any unitary 
matrix. Similar to coherent STBC case, when riR = 1, the op¬ 
timum power loading strategy is identical to the “water-filling” 
in information theory. When n u > 1, a generalized water- 
filling strategy gives the optimum P d . Following Sections give 
the optimum pi for d27b for n r = 1, 2, 3 receive antennas. For 
other cases, the the generalized method discussed in Section 
IIY-El can be directly applied to obtain the optimum pi for C3- 


B. MISO Channel 


The optimization problem involved in this case is similar to 
the water-filling problem in information theory, which has the 
optimal solution 


Pi = 


Vd <ti, 

0 , otherwise, 


(28) 


where the water-level l/v d is chosen to satisfy 



7 firiT 

8 + (3 


C. riT x 2 MIMO Channel 

The optimum pi for this case is 


Pi = 


A + \/~K, v d < U(ri + r 2 ); 
0 , otherwise. 


where v is chosen to satisfy 

tit 

max ^0, A + 

»=1 


"ffiriT 

8 + [3 


with 

A _ 2 rir 2 ff - v d tj(r 1 +r 2 ) 
2 v d r!r 2 ti 

and 


K = 


v d ti(ri 


- r 2 ) 2 + 
2v d rir 2 t\ 

























D. utxS MIMO Channel 

For the case of tit transmit antennas and n R = 3 receive 
antennas, the optimum pi is given by 


Pi = 


-j^ + Z, v d <ti(ri+r 2 +r 3 ); 

0, otherwise. 


where Vd is chosen to satisfy 



7 (3ht 

8 Tp' 


with 


Z = 
Zi = 


Z 2 + \J' zf + z\ 

3z!Z 3 - z% 

--> z 2 



9 , 272:223 — 27z 3 z\ — 2 
54z| 


z 3 = v d rir 2 r 3 t f, 2 2 = v d t\{r x r 2 + rir 3 + r 2 r 3 ) - 3rir 2 r 3 ff, 
21 = v d U(ri + r 2 + r 3 ) - 2tj (r 3 r 2 + rir 3 + r 2 r 3 ) and z 0 = 
Vd - U(ri + r 2 + r 3 ). 


E. Spatially Uncorrelated Receive Antennas 

If n rt receive antennas are placed ideally within the 
receiver region such that the spatial correlation between 
antenna elements is zero (i.e., J^ R Jr = I), then the objective 
function in reduces to a single summation and the 

optimum pi is given by the water-filling solution d} 
obtained for the MISO channel. This is not to say that such 
an ideal placement is possible even approximately. A similar 
result holds for the coherent STBC case. 


VI. Simulation Results: Coherent STBC 

In this section, we will illustrate the performance im¬ 
provements obtained from coherent STBC when the spatial 
precoder F c derived in Section TlV’-Al is used. In particular, the 
performance is evaluated for small antenna separations and 
different antenna geometries at the transmitter and receiver 
antenna arrays, assuming an isotropic scattering environment 
(independent and identically distributed entries in scattering 
channel matrix H 5 ). In our simulations we use the rate-1 
space-time modulated constellation constructed in [5] from 
orthogonal designs for two transmit antennas. Also use the 
rate 3/4 STBC code for tit = 3,4 transmit antennas given in 
[5]. Modulated symbols c(k) are drawn from the normalized 
QPSK alphabet {±l/y/2 ± i/y/2}. 

A. MISO Channels 

First we illustrate the water-filling concept for tit = 2,3 
and 4 transmit antennas, where the transmit antennas are 
placed in uniform circular array (UCA) and uniform linear 
array (ULA) configurations 6 with 0.2A minimum separation 
between two adjacent antenna elements. For each transmit 
antenna configuration we consider, Table-|J lists the radius 
of the transmit aperture, number of effective communication 

6 This precoder can be applied to any arbitrary antenna configuration. 


TABLE I 

Transmit antenna configuration details corresponding to 

WATER-FILLING SCENARIOS CONSIDERED IN FlG.0 


Antenna 

Configuration 

Tx aperture 
radius 

Num. of 
modes 

rank( J? J^ T ) 

2-Tx 

0.1 A 

3 

2 

3-Tx UCA 

0.115A 

3 

3 

3-Tx ULA 

0.2A 

5 

3 

4-Tx UCA 

0.142A 

5 

4 

4-Tx ULA 

0.3A 

7 

4 


modes 7 [18] in the transmit region and the rank of the transmit 
side spatial correlation matrix JtJ't- Note that, in all spatial 
scenarios, we ensure that JtJ't is full rank in order that the 
average PEP upper bound <1121 to hold. 
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Fig. 1. Water level {l/v c ) for various SNRs for a MISO system, (a) tit = 2, 
(b) n T = 3 - UCA, (c) n T = 4 - UCA, (d) n T = 3 - ULA and (e) n T = 4 
- ULA for 0.2A minimum separation between two adjacent transmit antennas. 


Fig. Qshows the water levels for various SNRs. For a given 
SNR, the optimal power value qi is the difference between 
water-level \/v c and base level 1 /t 7 , whenever the difference 
is positive; it is zero otherwise. Note that, with this spatial 
precoder, the diversity order of the system is determined by 
the number of non-zero <&’s. It is observed that at low SNRs, 
only one qi is non-zero for tit = 2 and 3-UCA cases. In these 
cases, all the available power is assigned to the highest eigen¬ 
mode of JtJ't (° r to the single dominant eigen-channel of 
H ) and the system is operating in eigen-beamforming mode. 
With other cases. Fig. Etc), (d) and (e), systems are operating 
in between eigen-beam forming and full diversity for small 
SNRs as well as moderate SNRs. In these cases, the spatial 
precoder assigns more power to the higher eigen-modes of 
JtJt (or to dominant eigen-channels of H) and less power 
to the weaker eigen-modes (or to less dominant eigen-channels 
of H). 

7 The set of modes form a basis of functions for representing a multipath 
wave field. 
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Fig. 2. Performance of spatial precoder with two transmit and one receive 
antennas for 0.2A separation between two transmit antenna elements: rate-1 
coherent STBC. 


Fig. 3. Performance of spatial precoder with three transmit and one 
receive antennas for 0.2A minimum separation between two adjacent transmit 
antennas for UCA and ULA antenna configurations: rate-3/4 coherent STBC. 


Fig.0 illustrates the BER performance of the rate 1 STBC 
with and without spatial precoder for nr = 2. It can be 
observed that at very low SNRs, we obtain a pre-coding gain of 
about 1.5dB. In fact, at very low SNRs, the optimum scheme 
is equivalent to eigen-beam forming. However, as the SNR 
increases, the precoder becomes redundant and the optimum 
scheme approaches STBC, where it operates in full diversity. 
This corroborates the claim that the 2x1 STBC has good 
resistance against the spatially correlated fading at high SNRs 
as shown in [24], 

BER performance results for 3-Tx UCA, ULA and 4-Tx 
UCA, ULA antenna configurations are shown in Fig. 0 and 0 
respectively for rate 3/4 STBCs. For 3-Tx UCA, the results 
obtained are similar to the results of nr = 2 case above. In this 
case, at low SNRs, the system operates in eigen beam-forming 
mode and at high SNRs, it is operating in full diversity mode as 
shown in Fig. 0b). For the other three cases, it is observed that 
the optimum scheme provides a clear performance advantage 
over the STBC only system for all SNRs concerned. For 
example, at 0.01 bit-error-rate, we obtain a precoding gain 
of about ldB. However, these systems operate in between 
eigen beam-forming and full diversity as the precoder assigns 
zero powers to some of the transmit diversity branches of the 
channel. As before, at higher SNRs, the system operates in 
full diversity and the optimum scheme approaches STBC. 

In all cases, at high SNRs we observed that ULA antenna 
configuration provides better performance than UCA antenna 
configuration when the spatial precoder is used. This is be¬ 
cause, the number of effective communication modes in the 
transmit region is higher for the ULA case (large aperture 
radius of ULA, c.f. Table Q} than the UCA case and the spatial 
precoder efficiently activates the modes in the transmit region 
of ULA. This observation suggests that our precoding scheme 
gives scope for improvement of ULA performance at high 
SNR. 



Fig. 4. Performance of spatial precoder with four transmit and one 
receive antennas for 0.2A minimum separation between two adjacent transmit 
antennas for UCA and ULA antenna configurations: rate-3/4 coherent STBC. 


B. MIMO Channels 

We now examine the performance of the spatial precoder for 
multiple transmit and multiple receive antennas. For example, 
we consider nr = 2,3 transmit antennas and ur = 2 
receive antennas. In all cases, two receiver antennas are placed 
A apart, which gives negligible effects on the performance 
due to antenna spacing. As before, the minimum separation 
between two adjacent transmit antennas is set to 0.2A. Note 
that this situation reasonably models the uplink of a mobile 
communication system. For each case, the optimum qi is 
calculated using (I19> . Fig- 0 illustrates the BER performance 
results for 2-transmit, 2-receive antennas for rate 1 STBC and 
Fig- |U illustrates the BER performance results for 3-transmit, 
2-receive antennas for rate 3/4 STBC. Performance results 
obtained here are similar to that of MISO cases above. 



















10° 


10° 


Coherent STBC 2 Tx 2-Rx 
O Coherent STBC with precoder: 2-Tx 2 Rx 



10 -5|-1-1-1- 1 - 1 - 1 - 1 - 

0 2 4 6 8 10 12 14 16 

SNR (dB) 

Fig. 5. Performance of spatial precoder with two transmit and one 
receive antennas: receive antenna separation A and minimum transmit antenna 
separation 0.2A for UCA and ULA antenna configurations: rate-3/4 coherent 
STBC. 



Fig. 6. Performance of spatial precoder with three transmit and two 
receive antennas: receive antenna separation A and minimum transmit antenna 
separation 0.2A for UCA and ULA antenna configurations: rate-3/4 coherent 
STBC. 


VII. Simulation Results: Differential STBC 


We now demonstrate the performance improvements ob¬ 
tained from differential space time block coded systems when 
the spatial precoder derived in Section IV-AI is applied. As 
before, the performance of differential space-time coded sys¬ 
tems is investigated for small antenna separations and different 
antenna geometries assuming a rich scattering environment 
surrounding the transmit and receive antenna arrays (i.e., i.i.d 
entries in Hs). We use the rate-1 space-time modulated 
constellations constructed in [5] from orthogonal designs for 
two and four transmit antennas. Normalized QPSK alphabet 
{±l/v/2± i/y/ 2} and normalized BPSK alphabet {±l/\/2} 
are used with two and four transmit antenna STBC, respec¬ 
tively. 


DSTBC scheme w o precoder 
* DSTBC scheme with precoder 
0 STBC coherent detection 



SNR (dB) 

Fig. 7. Performance of spatial precoder with two transmit and one receive 
antennas for 0.1 A separation between two transmit antennas: rate-1 differential 
STBC. 

A. MISO Channel 

Fig □ illustrates the BER performance of the differential 
STBC with and without spatial precoder when nr = 2. Also 
shown for comparison is the BER performance of the STBC 
when coherent detection is employed at the receiver. In all 
cases, two transmit antennas are placed 0.1 A distance apart. 
It can be seen that at the BER of 0.05, the performance of 
the precoded system is 1.25dB better than that of the non- 
precoded differential orthogonal space-time coded system and 
1.75dB away from the coherent detection case. However at 
high SNRs, the precoder becomes redundant and the optimum 
scheme approaches differential STBC. 

BER performance results for 4-Tx UCA and 4-Tx ULA 
antenna configurations are shown in Fig. |U and El respectively. 
For both antenna configurations, the minimum separation 
between two adjacent antenna elements is set to 0.2A, cor¬ 
responding to aperture radii 0.142A and 0.3A for UCA and 
ULA antenna configurations, respectively. Simulation results 
show that the BER performance of the optimum scheme is 
better than that of the differential STBC system for both 
antenna configurations. For example, at 10 2 BER, we obtain 
precoding gains of about ldB and 1.5dB with UCA and ULA 
antenna configurations, respectively. In comparison with the 
coherent detection at the receiver, BER performance of the 
optimum scheme is 2dB and 1.5dB away for UCA and ULA 
antenna configurations, respectively. 

B. MIMO Channel 

We now examine the performance of the proposed optimum 
scheme for multiple transmit and multiple receive antennas. 
As an example, we consider a MIMO system consisting of 
tit = 2 transmit antennas and n r = 2 receive antennas. 
The two receiver antennas are placed A apart, which gives 
minimum effect on the performance due to antenna spacing at 
the receiver antenna array, and the two transmit antennas are 
placed 0.1 A distance apart. Note that this situation reasonably 
models the uplink of a mobile communication system. Fig. 
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Fig. 8. Performance of spatial precoder with four transmit and one Fig. 10. Performance of spatial precoder with two transmit and two receive 
receive antennas for 0.2A minimum separation between two adjacent transmit antennas. Transmit antenna separation 0.1 A and receive antenna separation A: 
antennas; UCA transmit antenna configuration: rate-1 differential STBC. rate-1 differential STBC. 



tures. Using 0 , we define the modal correlation between 
complex scattering gains as 

Assume that the scattering from one direction is independent 
of that from another direction for both the receiver and the 
transmitter apertures. Then the second order statistics of the 
scattering gain function g(q i>, ip) can be defined as 

E {g(4>, <P)9*(<!>', <f')} = - <t>')S(ip - tp'), 

where G(cj>,ip) = E {\g{<j>, <^)| 2 } with normalization 

ff G{(j>,ip)d(pd(j) = 1. With the above assumption, the modal 

t i' 

correlation coefficient, 7 ^ m , can be simplified to 


Fig. 9. Performance of spatial precoder with four transmit and one 
receive antennas for 0.2A minimum separation between two adjacent transmit 
antennas; ULA transmit antenna configuration: rate-1 differential STBC. 


M shows the performance of the optimum scheme with 
two transmit and two receive antennas. Performance results 
obtained here are similar to that of MISO cases considered 
above. 

Note the objective function of D-STBC optimization prob¬ 
lem is derived for high SNR. However, from our simulation 
results, we observed that proposed precoding scheme also 
gives good performance at low SNRs. 


C. Effects of Non-isotropic Scattering 

In practise, wireless channels experience non-isotropic scat¬ 
tering (limited angular spread about a mean angle of depar¬ 
ture/arrival) both at the transmitter and the receiver antenna 
arrays. We now investigate the effects of non-isotropic scat¬ 
tering on the BER performance of differential STBC when the 
spatial precoding scheme derived in Section IV-AI is used. 

First we derive expressions for correlation between different 
communication modes at the transmitter and receiver aper¬ 


SixS 1 


ip)e 




^dipdtj). 


Then the correlation between f-th and f'-th modes at the 
receiver region due to the m-th mode at the transmitter region 
is given by 


l 1 ' 1 ' = [ V Rx (ip)e e ' )v d(p, V m, (29) 
J§ i 


where Vr x (<p) = fG((j>,ip)d(t> is the normalized azimuth 
power distribution (APD) of the scatterers surrounding the 
receiver antenna region. Here we see that modal correlation 
at the receiver is independent of the mode selected from 
transmitter region. 

Similarly, we can write the correlation between m-th and 
m-th modes at the transmitter region due to the /:-th mode at 
the receiver region as 


]m,m' 




V£, 


(30) 


where Vtx(4>) = fG((f>,<p)d(p is the normalized azimuth 
power distribution at the transmitter region. As for the receiver 
modal correlation, we can observe that modal correlation at the 
transmitter is independent of the mode selected from receiver 



















region. Note that, azimuth power distributions Vr x (tp) and 
V Tx {0) can be modeled using all common power distributions 
such as uniform-limited [25], Gaussian [26], Laplacian [27], 
cos 2p <f> distribution [25], etc. 

Denoting the p -th column of scattering matrix Hr as Hs, P , 
the (2 Nr + 1) x ( 2Nr + 1) receiver modal correlation matrix 
can be defined as 

Mr±e{h 3 , p hI p }, 

where (£,£')-th element of Mr is given by < 1291 above. 
Similarly, the transmitter modal correlation matrix can be 
defined as 

where Hs, q is the q-th row of Hr- (m,m')-th element of 
Mr is given by (l30l and M T is a (2 N T + 1) x (2 N T + 1) 
matrix. 

1) Kronecker Model as a Special Case: The correlation 
between two distinct modal pairs can be written as the product 
of corresponding modal correlation at the transmitter and the 
modal correlation at the receiver, i.e., 

7m,m' = 7 ’ 7 m,m'- (31) 

Facilitated by m. we write the covariance matrix of the 
scattering channel if s as the Kronecker product between the 
receiver modal correlation matrix and the transmitter modal 
correlation matrix, 

Rs = E 1 = Mr ® Mr . (32) 

Note that (1311 holds only for class of scattering environments 
where the power spectral density of modal correlation function 
satisfies [28,29] 

G{ct>^) =V Tx {ct>)V Rx {p). (33) 

Note that, 1331 is the necessary condition in which a channel 
must satisfy in order for (1321 to hold . 

Assuming Rr is a positive definite matrix, a channel 
realization of the scattering channel Hr can be generated by 

vec(H s ) = R 1 s /2 vec(W s ), (34) 

where R s ! is the positive definite matrix square root [22] 
of Rr and W r is a (2Nr + 1) x (2 Nr + 1) matrix which 
has zero-mean independent and identically distributed complex 
Gaussian random entries with unit variance. Furthermore, 
using J32> . the full correlation matrix of the MIMO channel 
H, given by 0, can be written as 

R = (j* R M R jty ® (j T M T 4) . (35) 

For simplicity, here we only consider the modal 
correlation at the transmitter region and assume the 
effective communication modes available at the receiver 
region are uncorrelated, i.e. Mr x = I- 2 N Ti +'\ ■ It was shown 
in [30] that all azimuth power distribution models give 
very similar correlation values for a given angular spread, 
especially for small antenna separations. Therefore, without 



Fig. 11. Precoder performance in non-isotropic scattering environments, 
at = 30° mean AOD </>o = 0° for a uniform-limited azimuth power 
distribution at the transmitter. 2x2 MIMO system. Transmit antenna 
separation 0.1 A and receive antenna separation A: rate-1 differential STBC. 


loss of generality, we restrict our investigation only to the 
uniform-limited azimuth power distribution, which is defined 
as follows: 


Uniform-limited Azimuth Power Distribution: When the 
energy is departing uniformly to a restricted range of az¬ 
imuth angles ±A around a mean angle of departure (AOD) 
(f >o £ [—7r, 7r), we have the uniform-limited azimuth power 
distribution [25] 


nt) 


i 

2A ’ 


l</>-<&>I < A, 


where A represents the non-isotropic parameter of the azimuth 
power distribution, which is related to the standard deviation 
of the distribution (angular spread <r t = A/\/3). For the above 
APD, the (m,m')-th entry of Mr is given by 

((m - m')A)e < <"‘-"*'>*>. 

Figures II 1 1 and [21 show the BER performance of rate-1 
differential STBC code with two transmit antennas for the 
spatial arrangement considered in Section lYII-BI for transmit¬ 
ter angular spreads a t = 30° and 10° about the mean AOD 
(f> o =0°. The channel is modeled using 0 and < 01 . 

From Figures^2 anc lElit is observed that in the presence of 
non-isotropic scattering at the transmitter, proposed precoding 
scheme provides significant BER improvements at low SNRs. 
To further improve the performance, following Section El 
a precoding scheme can be easily derived by including the 
non-isotropic scattering parameters (angular spreads and mean 
AOA/AOD) at both ends of the MIMO channel. Unlike in 
the fixed precoding scheme, modified scheme will require 
the receiver to estimate and feedback scattering distribution 
parameters to the transmitter whenever there is a change in 
these parameters. 


VIII. Performance in other Channel Models 

Simulation results presented in previous sections used the 
channel model H = JrHrJ' t , which is derived based 











Fig. 12. Precoder performance in non-isotropic scattering environments, 
at = 10° mean AOD <j >o = 0° for a uniform-limited azimuth power 
distribution at the transmitter. 2x2 MIMO system. Transmit antenna 
separation 0.1 A and receive antenna separation A: rate-1 differential STBC. 


on plane wave propagation theory, to simulate the underly¬ 
ing channels between transmit and receive antennas. In this 
section we analyze the performance of precoding schemes 
(coherent and differential) derived in this paper applied on 
other statistical channel models proposed in the literature. 
In particular we are interested on channel models that are 
consistent with wave propagation. MISO and MIMO channel 
models proposed by Chen et al. [31] and Abdi et al. [32], 
respectively are two such example channel models. Sections 
IVIII-AI and IVIII-BI provide simulation results of coherent 
STBC applied on Chen’s MISO channel model and differential 
STBC applied on Abdi’s MIMO channel model, respectively. 
In following simulations, precoders are derived using Jt 
and Jr for given antenna configurations and the underlying 
channel H is simulated using Chen et al. and Abdi et al. 
channel models. 


A. Chen et al.’s MISO Channel Model 

Fig - El depicts the MISO channel model proposed by Chen 
et al., where the space-time cross correlation between two 
antenna elements at the transmitter is given by 


[-RML, „ = exp 


j~r~(dm - dn) 
A 


(36) 


Jo 


2n ^J [for cosy + "xO + [for sin'r - 


with 

mn — ~t “ j \^ J rnn Jn ) COS O-mn COS f3rnn\ 5 

1 CLn 

zmn ~j 7 j {dm d n { cos cr mn sin/3 mrl , 

dm “r d n 

a is the scatterer ring radius, 7 is the moving direction of 
the receiver with respect to the end-fire of the antenna, J'd 
is the Doppler spread and d mn is the receiver distance to the 
center of the transmit antenna pair m, n. All other geometric 
parameters are defined as in Fig. El 



Fig. 13. Scattering channel model proposed by Chen et al. for three transmit 
and one receive antennas. 


Fie. 1 1 4l shows the performance of spatial precoder derived in 
Section ITV-AI for rate-3/4 coherent STBC with three transmit 
antennas placed in a ULA configuration. In this simulation, 
we assume the time-varying channels are undergone Rayleigh 
fading at the fading rate fr>T = 0.001, where T is the 
codeword period. We set parameters a = 30A, d\ 2 = 
d s 2 p 3 = 0.2A, di 2 = 1000A, 7 = 20° and /3 U2 = 60°. 
All other geometric parameters of the model in Fig. 1 1 31 can 
be easily determined from these parameters by using simple 
trigonometry. In this simulation, a realization of the underlying 
space-time MIMO channel is generated using 03 and 03- 
From Fig. El we observed that proposed spatial precoding 
scheme gives significant performance improvements for time- 
varying channels. For example, at 0.05 BER, performance of 
the spatially precoded system is ldB better than that of the 
non-precoded system. 


B. Abdi et al.’s MIMO Channel Model 

In this model, space-time cross correlation between two 
distinct antenna element pairs at the receiver and transmitter 
is given by 

exp[jc P g cos(apq)] w 

Iq ({ k —a — bi m — c p q A sin (cx P q) 

2, abim cos (film / y) sin(o;pg) 

x [osin( 7 ) - bi m sin (Pi 

m)\ 

- j2n [a cos(/r - 7 ) - h m cos (/r - # m ) 

- c pq Asm(a pq )sm(n) ])} 1/2 ) , ( 37 ) 

where a = 2itf D r, b im = 2ndi m /X , c pq = 2n6 pq /X; f D is the 
Doppler shift; // is the mean angle of arrival at the receiver; 
































Fig. 14. Spatial precoder performance with three transmit and one receive 
antennas for 0.2A minimum separation between two adjacent transmit an¬ 
tennas placed in a uniform linear array, using Chen et al’s channel model: 
rate-3/4 coherent STBC. 



Fig. 16. Spatial precoder performance with two transmit and two receive 
antennas using Abdi et al’s channel model: rate-1 differential STBC. 



Fig. 15. Scattering channel model proposed by Abdi et ah for two transmit 
and two receive antennas. 


k controls the spread of the AOA; and 7 is the direction of 
motion of the receiver. Other geometric parameters are defined 
in Fig. El Note that this model also captures the non-isotropic 
scattering at the transmitter via A and the model is valid only 
for small A [32], 

Fig- El shows the performance of spatial precoder derived 
in Section IV-AI for rate-1 differential STBC with two transmit 
and two receive antennas for a stationary receiver (i.e. fn = 
0). In this simulation we set S 12 = 0.1A, d \2 = A and 
a \2 = P12 = 0°. We assume the scattering environment 
surrounding the receiver antenna array is rich, i.e., k = 0 
and the non-isotropic factor A at the transmitter is 10°. We 
assume the scattering channel satisfies the power distribution 
condition ©. A realization of the underlying MIMO channel 
is generated using d34b and J37l> . It is observed that our 
precoding scheme based on antenna configuration details give 
promising improvements for low SNR when the underlying 
channel is modeled using Abdi’s channel model. Therefore, 
using the previous results from Chen’s channel model and 
the current results, we can come to the conclusion that our 
fixed spatial precoding scheme can be applied to any general 
wireless communication system. 


IX. Concluding Remarks 


In this paper, by exploiting the spatial dimension of a 
MIMO channel we have proposed spatial precoding schemes 
for coherent and differential space-time block coded systems. 
Precoders are derived by minimizing certain upper bounds for 
the PEP subject to a transmit power constraint and assuming an 
isotropic scattering environment surrounding the transmit and 
receive antenna arrays. The proposed precoders are designed 
based on previously unutilized fixed and known parameters of 
MIMO channels, the antenna spacing and antenna placement 
details. Therefore, with these schemes the transmitter does not 
require any feedback of channel state information from the 
receiver, which is an added advantage over the other precoding 
schemes found in the literature. Since the precoder is fixed 
for fixed antenna configurations, proposed precoding schemes 
can be applied in non-stationary scattering channels as well as 
stationary scattering channels. 

We showed that proposed precoding schemes reduce the 
detrimental effects of non-ideal antenna placement and im¬ 
prove the performance of space-time coded MIMO systems. 
Precoders achieve these performance improvements by vir¬ 
tually arranging antennas into an optimal configuration as 
such the spatial correlation between all antenna elements is 
minimum. For 1-D arrays (ULA), we observed that precoder 
gives scope for improvement at high SNRs, but for 2-D arrays 
(UCA), improvements are only seen at low SNRs. 

Although the proposed precoders are derived for isotropic 
scattering environments, we observed that these precoders give 
significant performance improvements in non-isotropic scat¬ 
tering environments. Based on the performance improvements 
we observed, we believe that proposed schemes can be applied 
on uplink transmission of a mobile communication system as 
the proposed schemes can effectively reduce the effects due 
to insufficient antenna spacing and antenna placement at the 
mobile unit. 















Appendix I 

Proof of PEP Upper bound: Coherent Receiver 

The conditional average pairwise error probability P(S',; —> 
Sj), defined as the probability that the receiver erroneously 
decides in favor of Sj when Si was actually transmitted for 
a given channel, is upper bounded by the Chernoff bound [3] 


P(Si 


'j\h) < exp ( 


Sj) 


(38) 


where d 2 (Si, Sy) = h[I 
SjKSj-SjFl h = (vec (H 1 )) 


S A ]h\ S A = F d (Si - 
h = (vec(ff r )) T a row vector and 
7 = E s /o 2 is the average SNR at each receiver antenna. 
To compute the average PEP, we average m over the joint 
distribution of h. Assume h is a proper complex 8 titUr- 
dimensional Gaussian random vector with mean 0 and co- 
variance matrix Rh = E j h ' h X , then the pdf of h is given 
by [33] 


p(h) = 


1 


■ exp { — hR^h^}, 


t^utur \R h \ 

provided that Rh is non-singular. Then the average PEP is 
bounded as follows 

1 


P (Si 


S ■) <- 

0) - n nTn R IR h \ 


exp{—hR Q 1 h ] }dh 


(39) 


where Rq 1 = (jI nR ® S a + R)h)- Assume Rh is non¬ 
singular (positive definite), therefore the inverse R r j is pos¬ 
itive definite, since the inverse matrix of a positive definite 
matrix is also positive definite [22, page 142]. Also note that 
S a is Hermitian and it has positive eigenvalues (through code 
construction, e.g. [3]), therefore S a is positive definite, hence 
I nR A S a is also positive definite. Therefore i?. () 1 is positive 
definite and hence Rq is non-singular. Using the normalization 
property of Gaussian pdf 

1 


7T n T n R |J? 0 | J 
we can simplify (1391 to 


exp{— hR 0 1 h) }dh = 1 


P(Si - Sj) < 

or equivalently 

P {Si - Sj) < 
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Rq 1 Rh I 
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Appendix II 

Proof of Generalized Water-filling Solution for 
nn = 2 Receiver Antennas 
Let n,R = 2 in m , then we obtain the second-order 
polynomial rir 2 v c t 2 q 2 + {v c U{r\ + r 2 ) - 2rir 2 t 2 )q l + 
(v c — r\ti — r 2 ti) in q which has roots i = A + \J~K and 
q l 2 = A — \f~K, where A and K are given by I20> . Then the 

8 To be proper complex, the mean of both the real and imaginary parts of 
Hs must be zero and also the cross-correlation between real and imaginary 
parts of Hs must be zero. 


product qi^q it2 = ( v c - nU - r 2 ti)lr\r 2 v c t\. 

Case 1: < 7 i,igi , 2 > 0 v c > U(r\ +r 2 ). In this case, both 
roots are either positive or negative. Let v c = ati(r i + r 2 ), 
where a > 1. Then A = — tfa[{ri + r 2 ) 2 — 2rir 2 /a] < 0 
for all oi 1. Since /( ^ 0, (/i j2 0, thus q %,l must also 
be negative to hold v c > ti(r i + r 2 ). Therefore, when 
Vc > +r 2 ), the optimum qi is zero to hold the inequality 

constraints of J16t . 


Case 2: qi A qi, 2 < 0 => v c < L(ri + r 2 ). In this case, we 
always have one positive root and one negative root. Assume 
qi t i > 0 and qi t2 < 0 and let v c = ati{r\ + r 2 ), where 
0 < a < 1. For g,; i to positive, we need to prove that \f~K > 
tfa[{r\ + r 2 ) 2 — 2rir 2 /a] for 0 < a < 1. Instead, we show 
that 


\f~K < tfa[(ri + r 2 ) 2 - 2rir 2 /a], (40) 

only when a > 1. Note that, since K > 0, (1401 can be squared 
without affecting to the inequality sign. Therefore squaring 
EB and further simplification to it yields a > 1. This proves 
that q, -[ > 0 and qi t2 < 0 when v c < U(ri + r 2 ) and the 
optimum solution to dl6b is given by q^i. 


Appendix III 

Proof of Generalized Water-filling Solution for 
nu = 3 Receiver Antennas 

Let tir = 3 in (II 8 bt . then we obtain the third-order 
polynomial a^qf + a 2 q? + a\qi + do in qi which has roots 
[34] 


Qi, i — ~ + (S + T), 

q i,2 = ~-\{ S + T )+ l -^{S- T ), 


where S ± T = 


r + VQ 3 + R 2 \ ± [R - VQ 3 + R2 

and all other variables are as defined in Section HWDl then 
the product qi,iqi, 2 qi ,3 = (nU + r 2 U + r 3 ti -v c )/rir 2 r 3 v c t?. 


Case 1: qi,iqi, 2 qi, 3 < 0 => v c > t t (r-i + r 2 + r 3 ). Let 
v c = atj (r'i + r 2 + r 3 ), where a > 1. For a > 1, it can be 
shown that ( Q 3 + R 2 ) > 0, hence < 0 and (?i, 2 ,<Z *,3 £ C. 
Therefore, when v c > ti(r\+r 2 +r 3 ), the optimum qi is zero. 


Case 2: q l ^q l . 2 q l . 3 > 0 => v c < U{r 1 + r 2 + r 3 ). Let 
v c = ati(ri + r 2 + r 3 ), where 0 < a < 1. For 0 < a < 1, 
it can be shown that ( Q 3 + R 2 ) < 0 and hence 

we get two negative roots (/, : 2 , q i 3 < 0 and one positive root 
<371 > 0 as the roots of cubic polynomial. Therefore, when 
v c < ti(r\ + r 2 + r 3 ), the optimum solution to (1161 is given 

by qi, i- 
















Appendix IV 

Proof of the conditional mean and the 

CONDITIONAL VARIANCE OF U = 2Re{w(k)A] j y jl (k - 1)} 
A. Proof of Conditional Mean 

Mean of u condition on the received signal y (k — 1) can 
be written as 

rn u \ y (k-i) = E |2Re {w(fc) Aj J y t (/c - 1) j | y(k - 1)|, 

= 2R e^E{w(k) \ y{k - 1)} A \ tj y\k - 1) j . 

(41) 


is the covariance of the noise vector n(k — 1 ) condition on 
y{k — 1). Using the minimum mean square error estimator 
results given in [35], we can write 


^- l n(k—l)\y(k—l) ^n(k— l),n(k— 1) 


- E 


t 


l(k— l).n(k— 1) y(k—l).y(k—l) y(k l),n(k 1); 


t _ 2 y' — 1 

n Zj y(k—l),y(k— 1) 


(49) 


Substituting for E y(k—i) y{k— i) in (|49) and then the result 
in (I48b> gives the conditional variance cr 2 , y ( k _iy 


Substituting w(k) = n(k) — n(k — 1)5,; and noting 
E {n{k) | y(k — 1)} = 0, (14 1 ) can be simplified to 

fttu\y(k— 1 ) 2 Re ^TYl n ^k—l)\y(k—l)^i( k 1 )) ’ 

= 2Re |m Tl (/._ 1 )|j /( fe_ 1) (J — 5,5])y t (fc - 1) j , 

(42) 

where fh^k-i^k-i) = E {n(k - 1) | y(k - 1)}. Using the 
minimum mean square error estimator results given in [35, 
Section 2.3], we obtain 

Wn(fc-i)|i/(fc-i) =E{n(k- 1)} + [y(k- 1 )-E{y(k- 1)}] 

X ^- , y(k-l),y(k-l)^- , V l A-l),n(k-l)i 

where 

s w(fc-i),»(fc-1) = E W( k - - !)} ) (43) 

= E s X(k - 1 ) ] R H X{k - 1) + o*I nTnR , 

and 

^- J y(k—l),n(k—l) = E ^y^ [k 1 )Tl[k 1 )} , 

= a n^n T n R - (44) 

Since E {n(k — 1)} = 0 and E {y(k — 1)} = 0, we have 
m n (k-\)\ y (k-i) = crl,y(k - 1) (45) 

x (E a X(k - 1 )'R h X{ k - 1 ) + o 2 n I ) 

Substituting 1451 for fh n fi--i)\y(k-i) i n <1421 gives the condi¬ 
tional mean fh u \ y (k- 1 )- 


Appendix V 

Proof of PEP Upper bound: Non-coherent Receiver 

At asymptotically high SNRs, the PEP condition on the 
received signal y(k — 1 ) is given by 

P (Si - S 3 | y(k - 1)) = Q (J^j ■ 

Now using the Chernoff bound 

QO) < ^ exp (-y) > 
the conditional PEP can be upper bounded by 

P(5, ; - S 3 | y(k - 1)) < \ exp j . (50) 

To compute the average PEP, we average over the joint 
distribution of y(k — 1). Assume y(k — 1) is a proper complex 
Gaussian random vector that has mean E {y(k — 1)} = 0 and 
covariance 

Ry(k-i) =E{y\k-l)y{k-l)}, 

= E s X(k-l) i R H X(k-l) + a 2 n I nTnR (51) 

If R y(k _!) is non-singular, then the pdf of y{k — 1 ) is given 
by 

p(y(k — 1)) = f2 y exp j-y(fc- 1 )R~^ k _ 1) y' ! (k - l)j, 


B. Proof of Conditional Variance 

Variance of u condition on the received signal y{k— 1) can 
be written as 

°l\y(k-l) = E {\\ U ~ ™u\y(k-l) \\ 2 \y{k~l)} (46) 

E ^ TTlu\y(k— 1)) '^ J u\y{k— 1)) | 2/(^ 1) J* 

After some straight forward manipulations we can show 

'll '^u\y(k— 1 ) 2Rc { (?7<(/c) 1) 7Tl ri ^k— 1 ) \y(k— 1 )] 

x 5 i )A]' ii y t (fc- 1 )}. (47) 

Substituting <|47> for u — TO u | y (fc_i) in d46t gives E3, 
shown at the top of the next page, where E n (k),n(k) = 

E {mf (fc)n(fc)} = a 2 1 and 

’En(k-l)\y(k-l) = 

E{\\n(k - 1 ) - fh^k-^iyik-i) II 2 ! y{k - 1)} 


where fl y = n nT ™ K / Averaging over the pdf 

of y(k — 1 ), we obtain 

P(5, - S 3 ) < 

y j exp {—y(k - l)Rf l y ] {k - 1)} d y(k - 1), 

(52) 

where 

R d 1 = R y{k-l) + g ^2 Di ’P 

Assume Rh is non-singular (positive definite). It can be 
shown that both R y (k-i) al 'd E>i.j are positive definite. There¬ 
fore, Rd is non-singular. Using the normalization property of 
Gaussian pdf 

| fld | / exp {-y(fc - 1 )Rf 1 y\k - 1)} d y(k - 1) = 1, 
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can simplify <1521 to 

P(SWS,)< 


Ifld 
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or equivalently 

P(Si - Sj) < 
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